Code and Data Repository for Skitka et al. (2022): Probing the Nonlinear Interactions of Supertidal Internal Waves using a High-Resolution Regional Ocean Model

source.tar: This is all of the python code that must be in the path for the python compute and plot scripts to run.
aux_data.tar: There are some grid files that need to be referenced by the postprocessor in this directory

Postprocessing directories:
all_boxes: This contains partially processed data averaged across the five subregions studied in the paper.
all_boxes_v1: This contains all_boxes data from the first version of this paper submitted to ArXiV/JPO in Nov. 2022.  It does not contain scale-separated definitions of ID or SP.  
boxes24: partially processed data averaged over boxes B and D
boxes135: partially processed data averaged over boxes A, C, and E

Intructions for all_baxes/res/*:
Configure mkscript.sh so that the hard-coded directories are correct and, when appropriate, aux_data refers to the contents of aux_data.tar/res.  Run mkscript.sh in each subdirectory of all_boxes, boxes24, or boxes135.  Then run the plot directories with python.  E.g. python plot_flux_decomp_m.py.  Figures should be saved into the "figures" directory in each resolution directory.

Different plot scripts available:
plot_flux_decomp_m.py: plots the integrated flux budgets shown in figure 5.  This data can also be printed out / saved to make bar plots like figure 7.
plot_bispectrum_m.py: plots all of the bispectra shown in the paper.
plotspectral_T_F_m_omega.py: plots all of the T(m,omega) plots shown in figure 10.

Instructions for all_boxes/Comparison
Make sure directories at the top of plot_* are correct
run python3 plot_*
